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Abstract 



We present the result of a systematic study of pileup phenomena seen in the X-ray Imaging 
Spectrometer, an X-ray CCD instrument, onboard the Suzaku observatory. Using a data set of observed 
P ' sources in a wide range of brightness and spectral hardness, we characterized the pileup fraction, spectral 

hardening, and grade migration as a function of observed count rate in a frame per pixel. Using the pileup 
fraction as a measure of the degree of pileup, we found that the relative spectral hardening (the hardness 
ratio normalized to the intrinsic spectral hardness), branching ratio of split events, and that of detached 
events increase monotonically as the pileup fraction increases, despite the variety of brightness and hard- 
ness of the sample sources. We derived the pileup fraction as a function of radius used for event extraction. 
Upon practical considerations, we found that events outside of the radius with a pileup fraction of 1% or 
3% are useful for spectral analysis. We present relevant figures, tables, and software for the convenience 
, of users who wish to apply our method for their data reduction of piled- up sources. 

Key words: instrumentation: detectors — X-rays: general — X-rays: individual (Cygnus X-1) 
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1. Introduction CCD pixel, or one of its immediate neighbors during the 

exposure time per frame. It causes a complicated and irre- 
Charge coupled devices (CCDs) have been playing an versible phenomenon called "pileup" . For example, when 
important role in the modern X-ray astronomy (Lumb et two photons with energies of Ei and E2 fall in the same 
al. 1991) since the SXT (Soft X-ray Telescope; Tsuneta pixel within an integration period, it is impossible to dis- 
^ . et al. 1991) mounted on YOHKOH. As a photon-counting tinguish them from one with an energy of Ei + E2. It af- 
detector, the first successful use of it in space was realized fects obtained spectra, images, and light curves. Incorrect 
with the SIS (Solid-state Imaging Spectrometer; Burke et evaluation of pileup can lead to wrong interpretations of 
al. 1991) on board the ASCA satellite. the data. This is especially the case for bright sources 

The fifth Japanese X-ray satellite Suzaku (Mitsuda et because their small statistical uncertainty can be ovcr- 
al. 2007) also carries the CCD devices, the X-ray Imaging whelmed by systematic uncertainties brought by pileup. 
Spectrometer (XIS; Koyama et al. 2007), located at the Therefore, it is important to understand pileup effects in 
foci of the X-ray Telescope (XRT; Serlemitsos et al. 2007) order to maximize the wide-band capability of Suzaku. 
modules, as well as a non-imaging hard X-ray instrument. In general, when pileup occurs, flux and spectra become 
the Hard X-ray Detector (HXD; Takahashi et al. 2007). lower and harder, respectively. The pileup effects of the 
The XIS covers soft X-rays from 0.2 to 12 keV, while the SISs of ASCA, have been studied by Ebisawa et al. (1996), 
HXD covers from 10 to 600 keV. The wide-band coverage Ueda et al. (1997), and Kotani et al. (2000). Recently, 
by a combination of these instruments is a key competence a detailed study for the CCD detector of Chandra, the 
of Suzaku. ACIS, as well as a generalized formula of the pileup, has 

Sources that have rich statistics in the HXD are usu- been compiled by Davis (2001). However, the telescopes 
ally bright in the XIS as well. In observations of such a used for the XIS and the ACIS are entirely different in the 
bright source, more than one photon can strike the same shape and size of the point spread function (PSF), which 
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adds non-linear complications to the effect. Therefore, 
it is difficult to simply apply methods developed for the 
ACIS to the XIS. Likewise, the way of estimating pilcup^ 
for the CCD detector of XMM- Newton, the EPIC, can not 
be directly applied to the XIS. Thus, studies on the pileup 
effects should be performed specifically for the combined 
system of XISs and XRTs. 

One of the practical solutions to mitigate the pileup in 
the XIS is the "core-exclusion" method, in which the cen- 
tral region of the PSF is excluded to avoid the CCD region 
severely affected by pileup. This method was adopted in 
the analysis of 4U 1630-472 (Kubota ct al. 2007), Cygnus 
X-1 (Makishima et al. 2008), GX 339-4 (Yamada et al. 
2009), and others. Based on these experiences, we released 
a recipe to apply the method for XIS data as an official 
document^. In the document, it is poorly stated how to 
estimate the size of the core to be excluded. This is often 
difficult to do by using data that are affected by pileup. 

In this paper, we present the results of a systematic 
analysis for several bright sources affected by pileup in 
the XIS. A brief summary on pileup characteristic for 
phenomenological specific to the XIS-I-XRT system is pre- 
sented in section 2. After describing target selection and 
data processing in section 3, systematic studies on pileup 
effects are presented in section 4. In section 5, we discuss 
a criterion that utilizes pileup fraction and show concrete 
examples for the application. The main results of the pa- 
per are summarized in section 6. 

2. Characteristic of Pileup Effects 

2.1. Elements that Determine Pileup Effects 

Pileup itself is intrinsically a simple phenomenon. 
However, it is complicated by the following realistic is- 
sues: 

• sizes of charge clouds and pixels, 

• grade selection both in orbit and on ground, 

• the PSF of the X-ray mirror, 

• the exposure time per frame, or "snap time" . 

Therefore, the pileup should not be separated from the 
instrumental issues, the X-ray mirror, and data-processing 
specifics. 

When an incident photon is photo-absorbed by a Si 
atom in the CCD pixel, photo electrons create a charge 
cloud that expands by diffusion along the way to the 
electrode, where the distribution of charges depends on 
both the absorption depth and the size of the pixel. For 
any CCD detectors, the distributed pattern of the charge 
cloud, "grade" , is used to distinguish between X-ray and 
non X-ray events. The pattern changes as pileup becomes 
more severe. As easily expected, single pixel events are 
less likely affected by pileup than multiple pixel events. 
If all charge clouds created by a single photon were much 

^ http:/ /xmm.esac.esa.int/sas/currcnt/documcntation/ 

threads / epatplot .shtml 
^ http: / / www. astro, isas .ac.jp/ suzaku / analysis /xis/ pileup / 

HowToCheekPilcup_vl.pdf 
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smaller than 1 pixel, then it would be much easier to re- 
construct events. However, in reality, when an incident 
photon falls into a border of pixels, the charge clouds ex- 
tend larger than 1 pixel, up to four pixels. Thus, we must 
consider issues specific to instruments to assess pileup ef- 
fects. 

2.2. Features of the XIS + XRT Systems 

In addition to detailed information on the XIS described 
in Koyama ct al. (2007), some features related to pileup 
effects in the XIS and XRT systems are summarized. The 
XIS system consists of four CCD detectors (XISO, 1, 2, 
and 3). XISO, 2, and 3 have front-illuminated (FI) CCDs, 
while XISl contains a back-illuminated (BI) CCD. The 
depletion layers of the FI and BI CCD are 65 fim and 
42 fj,m, respectively. XIS2 has been dysfunctional since 
November 2006. Among the remaining CCDs, XIS3 is al- 
most identical to XISO, while XISl is different. Therefore, 
we treat XIS0(3) and XISl differently throughout this pa- 
per, unless otherwise noted. 

Events with grade 0, 2, 3, 4, and 6 are regarded as X- 
ray events, while events with other grades as cosmic-ray 
events (see Koyama et al. 2007 and figure 4 of Yamaguchi 
et al. 2006). The grade 1 and 5 events are charge dis- 
tribution called "detached events". As pileup becomes 
more severe, more X-ray events are regarded as grade 1, 
5, and 7, because of the overlaid distributions of multi- 
ple charge clouds. Note that, in the case of the normal 
clocking mode, the XIS has three editing modes, 2x2, 3x3 
and 5x5, which means that the pulse height values of the 
number of pixels are sent as a telemetry. The event de- 
tection algorithm is, however, common to all the editing 
modes; the central pixel and the surrounding 8 pixels are 
used to detect an event on orbit. On the other hand, the 
grading algorithm utilizes not only the 3x3 pixels of the 
event but also outer pixels, whose details depend on the 
editing modes. Note that most X-ray events, which are 
photo-absorbed not in neutral region but in a depletion 
layer, are not split into a region larger than 2x2 pixels. 

The pileup effects primarily depend on one parameter, 
the number of counts per snap time in a detection cell. In 
the case of the XIS, the number of detection cells is nine 
(3x3); i.e., the central pixel and the surrounding 8 pixels 
are used to detect the event in orbit. Efforts to mitigate 
pileup in the XIS are thus divided largely into two: (1) 
to reduce the effective frame time (snap time) by partial 
read of the CCD and (2) to use events in the outskirts 
of the PSF, in which the surface brightness is lower than 
that at the core. 

Each XIS has two clocking modes: P-sum and Normal. 
Throughout this paper, we treat data taken only with the 
latter mode. Snap time is basically 8 s in the Normal 
clocking mode. The XIS has two options that make the 
snap time shorter than 8 s to reduce counts per a frame 
time. One of them is the window option, in which only a 
certain part of the CCD, called window, is used for read- 
out. The other is the burst option, in which the snap time 
can be shortened arbitrarily by discarding events arriving 
in a certain period during a frame time. An interval be- 
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tween the end of the snap time and the start of the next 
one is caUed "delayed time" . The delayed time is s when 
the burst option is off, while it becomes longer than s 
when the burst option is employed. The window and burst 
options can be applied independently for each sensor. 

The combination of the XIS and the XRT provides 
a spatial resolution of ~ 2' in a half-power diameter 
(HPD). Note that one pixel of the XIS corresponds to 
1".042. The field-of-view (FOV) is 18' x 18' at the max- 
imum, which becomes smaller when a window option 
is used; 18' x (1/4 x 18') for 1/4 window option, and 
18' X (1/8 X 18') for 1/8 window option. Due to the at- 
titude fluctuation and calibration uncertainties, the 1/4 
window is recommended over the 1/8 window option. In 
many cases, a combination of the 1/4 window and burst 
options is used to mitigate the pileup. A rule of thumb for 
the count rate to cause a possible pileup for point sources 
for the Normal full window mode is ^12 s~^ (see 7.9.1 of 
Suzaku technical description'^). 

The XIS-I-XRT system has two features that make 
the core-extraction method practically viable; (1) energy- 
independence of the PSF shape (Serlemitsos et al. 2007), 
and (2) heavy over-sampling. In other words, we can per- 
form experiments of changing the surface brightness just 
by employing the different extraction radius from the PSF 
core, keeping the other instrumental responses unchanged. 
The CCD detectors onboard Chandra and XMM-Newton 
lack either or both of the two features, which makes the 
application and validation of the method more compli- 
cated. 

2.3. Observed Features Caused by Pileup Effects 

As a result of convolving mirror and detector responses 
with piled-up events, the following observational features 
appear. 

Image 

The observed peak of the PSF is smeared out. The 
smearing is often difficult to distinguish from the PSF 
core flattening due to the thermal wobbling (Uchiyama et 
al. 2008). In some cases, a hole at the center is created 
by the two reasons: (1) an event is regarded as an non 
X-ray grade due to many events falling around it, and (2) 
when, pileup becomes more severe, a sum of pulse height 
of piled-up events exceeds an upper threshold used to de- 
termine an event, and it is not regarded as an event. 
Spectra 

The measured count rate decreases in lower energy 
ranges, while it increases in the higher energy ranges, 
which artificially creates a hard tail in the spectrum. Since 
this feature is convolved with the intrinsic spectra, it ap- 
pears more clearly in a source with softer spectra than one 
with harder spectra, 

Light curves 

When the flux of the source changes during an observa- 
tion, the light curve can become plateaucd at the peak. 



Grade 

The branching ratio among the grades can change. In 
general, pileup decreases single-pixel events (grade 0) and 
increase multiple-pixel events (grade 2, 3, 4, 6, and 7). 

The extents of these phenomena depend on incident flux 
and spectra of a source. Based on these observational 
clues, we construct reliable criteria to judge and exclude 
the pileup effects. 

2.4. Definition of the Pileup Fraction 

In order to characterize the degree of pileup, we use 
the pileup fraction /pi following Davis et al. (2001). The 
definition of /pi is the ratio of P{k > 2,x) to P{k > l,x), 
where P{k,x) is the Poisson distribution function, x is 
defined as the mean counts per snap time per detection 
cell, and k is the actual number of photons falling per snap 
time per detection cell. It is expressed as 



_ E'^^^P{k,x) _ 1 
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x^ + 0{x^). 



(1) 



We expanded the equation by the second order of x since 
X is in many cases smaller than one. 

The theoretical value of f^ can be derived by substi- 
tuting X with the incident count. However, when pileup 
occurs, the incident count is unknown from the data. 
Therefore, we use the observed piled-up counts (xq^s) and 
an observed pileup fraction f°^^ as a proxy for a theoret- 
ical one f^. 

3. Data Set and Processing 

3.1. Target Selection 

To perform a systematic study on pileup, we need to se- 
lect as many point sources as possible to cover wide ranges 
of brightness of photons and intrinsic hardness, because 
the extent of pileup is determined by the former, while 
the appearance of pileup effects is sensitive to the latter. 

The list of the targets, as well as observational param- 
eters of these sources (e.g., the snap and delayed times) 
are summarized in tabic 1 . These meet the following crite- 
ria: (1) brighter than ^ 1 mCrab and (2) small variability 
with a root mean square of the count rate smaller than 
30% of its mean. The latter criterion is needed to avoid 
highly varying sources, in which pileup effects are difficult 
to disentangle with varying flux and spectra. 

3.2. Data Processing and Event Screening 

We performed data processing from unscreened event 
flies, and flrst discarded the telemetry-saturated periods 
by using the house keeping data* provided by the XIS 
team. The telemetry saturation itself does not change the 
spectral shape, though it may lead to underestimation of 
flux, and hence causing wrong estimation of pileup in some 
cases. 

We further screened the events in the same way as the 
standard criteria based on e.g., conditions on the satellite 



http: / /www. astro. isas.ac.jp/suzaku/doc/suzaku_td/ 
nodel0.html#SECTION001091000000000000000 



The file aeNNNNNNNNNxi [0-3] _0_tel_uf . gtl can be found in the 
distributed processing products. 
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Table 1. The basic information and the operated modes of the XISs for the selected samples. 



Nanie^ 


U DSld 


Date 




L (ks) 


XISO 




XISl 


XIS 3 














xicLte 


X 


b/U (sj 








GX 17+2 


402050020 


2007-09- 


27 


56.3 


0.5/1.5 


327 


5 


2.0/6.0 


346.6 


0.5/1.5 


375.5 


Serpens X-1 


401048010 


2006-10- 


24 


81.3 


1.0/1.0 


259 





1.0/1.0 


281.4 


1.0/1.0 


263.1 


f~\ -\r o Or\ A 

GX 339—4 


401068010 


2007-02- 


12 


256.5 


0.3/1.7 


950 


5 


0.3/1.7 


936.3 


0.5/1.5 


866.7 


Grab 


101004020 


2006-04- 


04 


40.4 


0.1/1.9 


1739 


6 


0.1/7.9 


2801.6 


0.1/1.9 


1710.6 


Cygnus X-l(#l) 


404075120 


2009-06- 


04 


44.5 


0.5/1.5 


412 


8 


0.5/1.5 


500.8 


0.5/1.5 


472.4 


ATT 1 f A A 

4U 1705—44 


401046020 


c\r\r\f^ r\c\ 

2006-09- 


18 


29.8 


2.0/0.0 


88 


4 


2.0/0.0 


86.2 


2.0/0.0 


91.0 


GX 349+2 


400003010 


2006-03- 


14 


58.4 


0.3/0.7 


339 


9 


0.3/1.7 


498.9 


0.3/0.7 


401.5 


Gygnus X-1 (#2) 


404075130 


2009-10- 


21 


/I A A 

40.0 


0.5/1.5 


228 


1 


0.5/1.5 


267.9 


0.5/1.5 


262.2 


4U 1630—472 


400010010 


2006-02- 


08 


43.8 


1.0/1.0 


112 





1.0/1.0 


116.2 


1.0/1.0 


127.8 


Anl X 1 




9007 00 


Zo 


1^ 7 


0.5/1.5 


159 


o 


0.5/1.5 


187.2 


0.5/1.5 


162.9 


Vela X-1 


403045010 


2008-06- 


17 


147.4 


2.0/0.0 


39 





2.0/0.0 


38.9 


2.0/0.0 


40.0 


3C 273 


702070010 


2007-06- 


30 


109.2 


8.0/0.0 


5 


1 


8.0/0.0 


7.6 


8.0/0.0 


5.5 


4U 1636-536 


401050010 


2007-02- 


09 


52.1 


1.0/1.0 


49 


2 


1.0/1.0 


56.7 


1.0/1.0 


48.7 


MCG -6-30-15 


700007010 


2006-01- 


09 


359.9 


8.0/0.0 


3 


9 


8.0/0.0 


6.6 


8.0/0.0 


3.7 


SS Cyg 


400006010 


2005-11- 


02 


80.8 


8.0/0.0 


3 


8 


8.0/0.0 


6.4 


8.0/0.0 


3.7 


MCG -5-23-16 


700002010 


2005-12- 


07 


219.6 


8.0/0.0 


3 


3 


8.0/0.0 


4.8 


8.0/0.0 


3.8 


PKS 2155-304 


104004010 


2009-05- 


27 


157.9 


2.0/0.0 


4 


8 


2.0/0.0 


9.7 


2.0/0.0 


5.1 



* Time from start to end of the observation. 

^ 'S and 'D' refer to snap time and delayed time in units of seconds, respectively. 
^ The count rate integrating over a whole region per seconds. 
Names are sorted by ascending order of counts per snap time in the XISl. 



orbit and data rates^. For the analysis of grade-branching 
ratio, we have extracted grade 02346, grade 0, grade 2346, 
and grade 1, and grade 5 events, separately. 
3.2.1. Attitude correction 
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Fig. 1. Time profiles of the image center in the SKY X and 
Y coordinate in the unit of pixel, taken from observation of 
Cygnus X-1 performed on October 21 in 2009. The bin size 
is 100 s. The open and filled circles are the values of SKY Y 
and X before and after applying the attitude correction with 
aeatttuned.py, respectively. 

As mentioned in Uchiyama et al. (2008), the attitude 
of the Suzaku satellite can be restored by using house 
keeping data including the temperature which has been 
implemented in the aeatt coord tool. Figure 1 shows the 
time profiles (open circles) of the barycenter position of 



Cygnus X-1 (#2) calculated from the observed image in 
the SKY X and Y (the coordinate corresponding to the 
position on the celestial sphere). Although the attitude 
was corrected with aeattcoord to some extent, it still 
shows a fluctuation of ^ 20 pixel in a time scale of 1 
ks. This fluctuation is non-negligible for evaluating the 
pileup effects at the peak of the PSF. 

We further corrected the image for the residual fluctu- 
ation using our python script aeatttuned.py, which has 
the same function as aeattcor2^ in the HEAsoft pack- 
age provided by John E. Davis, M. Nowak, and others in 
MIT. This method (1) creates the time profile of the im- 
age center of SKY X and Y in a time bin of 100 s, (2) 
modifies the attitude file to correct for the fluctuation, 
and (3) reprocesses the event file with the modified atti- 
tude file using xiscoord. We show the time profiles (filled 
circles) in figure 1 obtained after applying this method. 
Through this process, the fluctuation of SKY X and Y 
decreases to 3 pixel. In the following analysis, we have 
applied the attitude correction for all the events unless 
otherwise noted. Note that the method functions prop- 
erly only when a count rate is high enough to track image 
centers for every 100 s. 
3.2.2. Background subtraction 

We used the observation of the North Ecliptic Pole 
(NEP) to estimate the background count rate, taken on 
October 16 in 2007 (OBSID 100018010), considering the 
difference in a snap time. The counts per snap time (8 s in 



See http:/ /www. astro. isas.jaxa.jp/suzaku/process/ 
v2changes / criteria_xis.html. 



http: / /heasarc. gsfc.nasa.gov/ftools/caldb/help/aeattcor2. html 
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the case of NEP) per pixel of XISO and XISl are 1.6 x 10"^ 
and 3.8 x 10~^ counts frame^^ pixel^^, respectively. This 
is much smaller than most sources in our targets, but can 
be non-negligible of ^-^10% of the surface brightness at the 
PSF outskirt of the dimmest sources in our list. 

Since the NXB of the XIS changes by a factor of 2 
(Tawa et al. 2008), the systematic error in the subtraction 
of the NXB is expected to be ^^20% at most. In the fol- 
lowing analysis, we have estimated backgrounds from the 
NEP for corresponding grade selections, otherwise stated. 
When creating a radial profile of the surface brightness, 
we subtracted values in the NEP from those of the sources 
in each radial bin. 

The Cosmic X-ray Background (CXB) was not sub- 
tracted from the spectra, since it is almost negligible of less 
than 0.1% of all the source signals in any energy ranges. 
Hot pixels are excluded from all of the events throughout 
the paper. 

4. Observed EfTects on Pileup 

Spectral Hardening and Grade Migration 

We first demonstrate two effects of pileup; i.e., spectral 
hardening and grade migration, using real data. We use 
a Seyfert galaxy MCG -6-30-15 and a Galactic black hole 
binary GX 339-4 as a representative case, respectively of 
little and significant pileup. 
J^.l.l. Spectral hardening 

As explained in subsection 2.2, the PSF is independent 
of the incident energy for the entire CCD field of view 
without pileup. As an example, we analyzed the spectra 
of MCG -6-30-15. Figure 2 shows ratios of XISO spectra 
extracted from several inner regions of the PSF against 
those from an outer region in the outskirts of the PSF, 
where pileup is negligible. Figure 2 (a) and (b) are for 
MCG -6-30-15 and GX339-4, respectively. The ratios are 
normalized to be 1.0 at 1 keV. 

We confirmed in figure 2 (a) that the spectral shapes 
between the 1' circle and 1-2' annulus regions agree with 
each other within ^ 5% in the flux in all energy bands. 
On the contrary, when spectra are affected by pileup as 
shown in figure 2 (b), the ratio between inner and outer 
regions becomes energy- dependent. 

It is clear how significantly the spectra differ with each 
other depending on the choice of the region. The pileup 
effect on spectral distortion can be noticed in the following 
two features: 1) moderate hardening over the whole en- 
ergy range, and 2) rapid hardening above ^ 7 keV, proba- 
bly due to sharp decrease in the quantum efficiency of the 
CCD. In addition to these two features, line-like features 
are noticed in the spectral ratio around the instrumental 
edges (Si K at 1.8 keV in the XIS and Au M at 2.3 keV in 
the XRT). We consider this feature is caused by either the 
ratio between the spectrum with smeared-edge features 
owing to pileup and the one without it or the change in 
the charge transfer efficiency under a high count rate. The 
latter is not a direct outcome of pileup, because, when a 
count rate is high, a part of the events works as sacrifi- 
cial charges to improve the charge transfer efficiency for 



others (Gendreau 1995). This shifts the edge structure in 
the energy spectrum to higher energies, which mimics a 
line-like feature when a spectral ratio is calculated. 
4- 1-2. Grade migration 

To confirm the changes in grade branching ratio caused 
by pileup, we compared grade-sorted spectra and grade 
branching ratios of GX 339-4 and those of MCG -6-30-15. 
We extracted five spectra consisting of grade 02346, grade 
0, grade 2346, grade 1, and grade 5, which we present for 
XISO and XISl in left and right panels in figure 3, respec- 
tively, for MCG -6-30-15 (upper panels) and GX 339-4 
(lower panels). To further clarify the energy dependence 
in grade branching ratio, we plot the spectral ratio against 
X-ray event grades (grade 02346). 

In the MCG -6-30-15 data, grade events (red) and 
grade 2346 events (cyan) are dominant over the entire en- 
ergy range, while the grade 1 (green) and grade 5 (blue) 
events are less than 1% for both sensors. The ratios be- 
tween grade and grade 2346 of XISO are different be- 
tween XISO and XISl. In XISO, a crossover is seen in the 
contribution of grade and grade 2346 events; grade is 
the largest below ~ 6 keV and grade 2346 above ^ 6 keV. 
In XISl, on the other hand, no such crossover is seen and 
the grade 2346 events are the largest reaching at about 
60-80% in all energy bands. The difference is qualita- 
tively understood by the fact that the distance between 
the point of photoelectric absorption in the depletion layer 
and the electrodes is longer for XISl than XISO for their 
opposite directions of illumination. 

The apparent differences between MCG -6-30-15 and 
GX 339—4 include (1) the crossover energy of grade and 
grade 2346 in XISO moves from 6 keV to 5 keV, (2) the 
ratio for grade 1 becomes larger in GX 339-4 (~5%), than 
that in MCG -6-30-15 (^1%); the trend can be seen for 
both XISO and XISl, and (3) the fraction of the grade 5 
event is larger in GX 339-4 than MCG -6-30-15, both for 
XISO and XISl. 

4-2. Radial Profiles 

After applying the attitude correction (section 3.2.1), 
we extracted cleaned events (grade 02346) from the entire 
CCD regions. We use the energy ranges from 0.5 to 10.0 
keV to maximize the photon statistics. Figure 4 shows 
radial profiles of surface brightness for XISO and XISl 
separately for each source. 

GX 17-f2 has the largest counts per snap time in XISl. 
It suffers a significant loss of counts at the PSF center, 
particularly in XISl, which results in a broader profile 
than the others. 

Multiplying the values of the vertical axis in figure 4a 
and b by 9(= 3 x 3), and then inserting them into x in 
equation (1), we obtained f^i", which is plotted in figure 
4 (c) and (d) for XISO and XISl. As long as the sur- 
face brightness is small, the pileup fraction is proportional 
to the surface brightness as expected from equation (1), 
Therefore, the profiles in figure 4 (a) and (b) are similar 
to those of (c) and (d), respectively. 
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Fig. 2. Ratio of XISO spectra extracted from several inner regions of the PSF against those from an outer region in the outskirts 
of the PSF, where pileup is negligible, (a) MCG -6-30-15. The inner region was taken from r = 0—60 pixels, while the outer region 
was from r = 60-120 pixels, (b) GX339-4. The inner regions were taken from r = 0-25, 25-60, 60-75, 75-90, and 90-120 pixels, 
while the outer region was from 120—240 pixels. 



MCG -6-30-15 Grade02346 GradeO Grade 2346 Grade Grade 5 



100 

XISO 

1 



> 

<u 

^ 10"^ 
fl 

3 

U 1" 

1 

.2 

0.01 
0.001 



XISl 



■TT+riTTf ,,,i.,.,r*~ . , , , I . , , , I . . . . I . , , . I ...,i,,.,i,,,!L.|il, 



■ftf 



0.7 1 2 
GX 339-4 



5 7 10 0.7 1 



5 7 10 




- -XISl 



*++ +V+i+++r if (PI, 



5 7 10 0.7 1 2 
Energy (keV) 



5 7 10 



Fig. 3. The detecter-response-convolved spectra and spectral ratio to the standard event (grade 02346) of MCG —6-30-15 (top) 
and GX 339—4 (bottom). The spectra of XISO and XISl are shown in left and right, respectively. The black, red, cyan, green, and 
blue spectra refer to grade 02346, grade 0, grade 2346, grade 1, and grade 5 events, respectively. 
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Fig. 4. (a) The radial profile of XISO of counts rate in a frame and a pixel for several sources; the relation between target names 
listed in table 1 and the styles of the mark arc shown below. Their backgrounds are subtracted in a matter mentioned in section 
3.2.2. The profiles are averaged cvcnly-spaccd in a linear scale in r < 100 pixel, and in a log scale in r > 100 pixel. Although the 
error bars in Y-axis are actually plotted, they arc smaller than the size of their marks, (b) The same figure as (a) except for XISl. 
(c) The same figure as (a) except for a function of the measured pileup fraction, (d) The same figure as (b) except for a function of 
the measured pileup fraction. The dotted lines and dashed lines refer to 3% and 1% of the measured pileup fraction, respectively. 



4-3. Pileup Effects Measured by Pileup Fraction 

The pileup fraction is the fundamental measure of 
pileup, while we cannot tell f^^ from /"f^ as mentioned in 
section 2.4. A more practical approach is to use the effects 
of pileup, such as spectral hardening and event migration 
as discussed in section 4.1, as a measure of pileup. Using 
the set of sources, we derive the quantitative relation be- 
tween these effects as a function of pileup fraction. 
4-3.1. Spectral hardening 

As shown in figure 2, we have confirmed that the pileup 
causes the spectral hardening. To study the spectral hard- 
ening more systematically, we calculated hardness ratio as 
the count rate ratio between 7.0-10.0 and 0.5-10.0 keV 
at different radii for all sources (table 1), and plotted it 
against f°l'^ in figure 5 (a) and (b) showing that our sam- 
ples include sources with various intrinsic hardness. As 
expected, the hardness ratio increases monotonically as a 
function of increasing pileup fraction. 

In order to compensate for the intrinsic spectral hard- 
ness of the sources, we normalized the hardness ratio at 
f°^^ of 0.1%, which is shown in figure 5c and d. Despite 



the wide range of intrinsic spectral hardness, the relative 
spectral hardening shows the same trend. The hardening 
becomes significant when f°^^ reaches ^ 1%, and increases 
by a factor of - 2 from Z"*^" of 0.1% to 10%. The degree 
of hardening is larger for XISO than XISl because XISO 
has a larger effective area in the higher energy band. 
4.3.2. Grade migration — split events 

Since the grade branching ratio is sensitive to pileup, it 
can also be a good criterion on pileup, as shown in figure 
3. In figure 5 (e) and (f), we compare with the "split 
event ratios", the ratio of the split events (grade 2346) 
to the standard events (grade 02346). The energy range 
from 7 to 10 keV is used because the split events are more 
imminent in the hard band (figure 3). 

As shown in 5e and f, the split event ratios in all the 
17 observations agree with each other within ~ 10%. The 
branching ratios are fairly constant below ~ 1% of f^^- 
As f°^^ becomes larger than ~ 1%, the split event ratios of 
XISO and XISl start to increase. Similarly to the changes 
in spectral hardness, the trend of the change of the split 
event ratio as a function of pileup fraction does not depend 




Fig. 5. (a) The profiles of tfie iiardness ratios between 7.0-10.0 keV and 0.5-10.0 keV of tiie XISO, as a function of pileup fraction. 
Tlie symbols are the same as used in figure 4. (b) the same as (a), except for XISl. (c) the same as (a), except for being normalized 
to 1.0 at a pileup fraction of 0.1 %. (d) the same as (c), except for XISl. (e) the ratio of split events (grade 2346) to the standard 
events (grade 02346), as a function of pileup fraction, (f) the same as (e), except for the XISl. (g) the same as (e), except for 
detached events (grade 1). (h) the same as (g), except for the XISl. 



No. ] 

on intrinsic hardness. The onset of the pileup effect can 

be seen at several % of f°f^. 

4-3.3. Grade migration — detached events 

Finally, we show in figure 5g and h the ratio between the 
detached events (grade 1) and the standard events (grade 
02346). Since, as shown in figure 3, grade 1 events are not 
significantly dependent on energy, we adopted the energy 
ranges of 3-10 keV. 

The profiles of XISO agree with each other within ^ 
10%, which is reasonably understood by the idea that 
pileup causes the gradual increase of the detached events. 
Meanwhile, the profiles of XISl in 3C 273 and 4U 1705-44 
deviate from the others by a factor of ~ 5. We examined 
the images of the two observations, and found that it was 
caused by an asymmetric feature in the image of grade 
1 events for some reasons. Aside from these outliers, the 
ratio of grade 1 events monotonically increases by orders 
of magnitude with the pileup fraction, in contrast to the 
profiles of relative hardness ratio and the split event ratio. 
The latter two appear to have a threshold in the pileup 
fraction to cause any changes (figure 5 c-f), whereas the 
detached event ratio appears to change from the lowest 
level of pileup (figure 5 g and h). Thus, the fraction of 
the detached events is most sensitive to the pileup frac- 
tion. 

5. Practical case study 

5.1. Pileup Fraction Criteria 

We have estimated the radius corresponding to 1% and 
3% of /pf"* (hereafter Ri% and i?3%, respectively) by using 
the data points in figure 4 (c) and (d), and summarized 
them in table 2. We also derived these values (in paren- 
theses in table 2) for the images before attitude correction. 
The estimate of the radii does not change so much, ~ 1% 
at most, whether the attitude correction discussed in sec- 
tion 3.2.1 was performed. The method is thus useful for 
images without the correction. 

We mask events within a certain radius, where the 
pileup effect is non-negligible. We use the radius for 1% 
and 3% pileup fraction (i?i% and i?3%, respectively) as 
representative cases of conservative and liberal choices. 
GX 339-4, for example, has i?3% and Ri% of 58.5 and 
101.4 pixels in the XISO, according to table 2. As shown 
in figure 2, the spectral ratio of 90-120 pixel (^ Ri%) 
to 120-240 pixel seems almost fiat. Thus, Ri% can be a 
reasonable value for judging the onset of the pileup. In 
general, the appropriate mask radius should be chosen 
depending on science. At a larger R, the pileup effects 
decrease. However, attention should be paid for system- 
atic uncertainties besides pileup effects. We see such an 
example in section 5.3. 

Figure 6 shows the relation between the measured 
counts per snap time with the values of i?3% and Ri% 
listed in table 2. The values of Ri% in XISO reaches al- 
most zero at ^ 10 counts per snap time. To make it pos- 
sible for general users to estimate the pileup extent before 
performing an observation, we quantified f°^^ as a func- 
tion of the measured counts per snap time, by fitting these 
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Table 3. The fitting results of and Ri% as a function 

of measured surface density. 



XISO XISl 





(i?3%)* 


iRi%r 


{R3%r 


{Ri%r 


at 


60.7±2.2 


107.3 ±3.8 


58.7±2.0 


121.0±4.9 


hi 


-10" 


-9.2 ±0.7 


-10" 


-7.5±0.7 


c§ 


45.7±5.9 


10" 


74.9 ±5.3 


10" 


X(d.o.f) 


10.8(8) 


6.4(12) 


14.3(7) 


12.9(8) 



* The fitting results of i?3% and respectively, 
t In an unit of pixel. 

i In an unit of 10""^ counts"-'^ frame pixel. 
^ In an unit of counts frame" ^pixel"^. 
" The values are unconstrained and are fixed to the 
value at the lowest bound. 

values with the following phenomenological equation: 

Radius(a;) = a(l — exp(6(a; — c))) [pixel], (2) 

where x is the measured counts per snap time, while a, 6, 
and c are free parameters. We took into account a system- 
atic error of 10% for the count rates, and ±5 pixel as a sys- 
tematic error for the radius, considering the intrinsic vari- 
ation of the sources and individual spectral shapes. The 
data that show the radius less than 10 pixel and measured 
count rate larger than 320 counts snaptime"^ are not used 
for fitting. In addition, only for the data with a pileup 
fraction of 3%, those that show a measured count rate 
smaller than 70 counts snaptime"^ are not used. Since 
the observed PSF of the Crab Nebula is slightly wider 
due to its angular size of ~ 1', its values are not used 
for fitting. All the fits became acceptable, and the resul- 
tant parameters are summarized in table 3. We plotted 
the ratios to the best-fit model in the bottom of figure 6, 
showing an agreement to the data within ~ 20%. 

A quick way to find suitable observation modes for a 
bright source, is to adjust window modes or burst options 
to realize that counts per snap time over 0.5-10 keV be- 
comes below 10 counts per snap time. After observing 
it, then you can approximately know the extent of the 
pileup by comparing the measured counts per snap time 
with the radius in figure 6, or equation (2). 

5.2. Application to a severely piled-up source 

To study uncertainties by using a region outside R^% 
and we have compared spectra of the XISO, taken 

from the observation of Cyg X-1 on June 2009 in table 
1, which is the most severely affected by pileup among 
the 25 observations of Cyg X-1 in the Low/Hard state 
(for details see Yamada 2011). As shown in table 2, the 
values of R-^% and Ri% in the XISO are 46 and 88 pix- 
els, respectively. Figure 7 shows the detector-response- 
convolved/deconvolved spectra accumulated from a whole 
CCD region, outside i?3% and outside 
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Table 2. The values of Ri% and -R3% after and before the attitude correction. 



Name 


XISO 




XISl 




XIS 3 








o /o 






J- /o 


GX 17+2 


39.0(38.9)T 


78.3(78.0) 


94.3(94.4) 


151.2(151.1) 


42.3(42.1) 


81.7(80.9) 


Serpens X-1 


55.5(55.4) 


97.9(97.9) 


58.5(57.7) 


106.0(106.6) 


54.8(54.8] 


99.6(99.6) 


GX 339-4 


58.4(58.5) 


101.6(101.4) 


53.6(53.7) 


110.8(111.0) 


76.4(76.5) 


117.3(117.2) 


Crab 


43.2(43.4) 


83.4(83.3) 


42.0(42.0) 


97.0(97.1) 


39.7(39.9) 


82.0(81.8) 


Cygnus X-l(#l) 


46.0(44.9) 


88.6(88.8) 


49.4(49.5) 


101.1(101.0) 


54.5(54.7) 


92.5(92.3) 


4U 1705-44 


41.2(41.1) 


80.2(80.4) 


35.3(35.2) 


84.9(83.5) 


40.5(40.2) 


80.3(79.9) 


GX 349+2 


34.3(34.2) 


61.2(61.3) 


28.7(28.6) 


77.4(77.4) 


36.6(36.5) 


63.4(63.5) 


Cygnus X-1 (#2) 


28.9(29.0) 


64.9(64.8) 


27.4(27.4) 


74.6(74.9) 


32.0(32.1) 


70.7(70.6) 


4U 1630-472 


28.7(28.7) 


65.1(65.0) 


20.5(20.4) 


63.5(63.4) 


30.4(30.2) 


69.2(69.4) 




16.3(16.6) 


50.9(51.1) 


1 1 fin 1 Q^ 

ii.Ol^i i.o 


00.0(^04.0 ) 


iO.ol^io. i 


oi. i (^ou.y ) 


Vela X-1 


8.4(6.4) 


53.0(53.1) 


-(-) 


40.2(39.3) 


7.2(3.8) 


53.4(53.7) 


3C 273 


5.0(7.2) 


27.5(29.0) 


19.1(20.5) 


27.3(28.2) 


3.4(5.4) 


28.9(29.1) 


4U 1636-536 


2.8(3.5) 


37.2(36.8) 


4.0(4.3) 


34.8(34.6) 


2.3(1.9) 


36.1(36.2) 


MCG -6-30-15 


-(-) 


19.3(19.4) 


-(-) 


19.1(19.5) 


-(-) 


14.6(15.0) 


SS CYG 


-(-) 


19.8(20.1) 


-(-) 


19.7(20.6) 


-(-) 


15.3(15.3) 


MCG -^5-23-16 


-(-) 


15.4(15.3) 


-(-) 


8.3(7.9) 


-(-) 


18.3(18.4) 


PKS 2155-304 


-(-) 


-(-) 


-(-) 


5.5(6.2) 


-(-) 


-(-) 



* The radius at 3% and 1% pileup fraction, in an unit of pixel. 

^ The values in parenthesis are obtained without the attitude correction as mention in subsection 3.2.1. 
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Fig. 6. (top) The relation between measured counts per snap time (products of snap times by count rates as listed in table 1) over 
0.5-10.0 keV and and K\%. The black and red points correspond to K\% and Kz% of XISO, respectively; v^hereas the green 

and blue ones correspond to those of XISl, respectively. The best-fit models are overlaid as dotted curves in the same color as the 
data. The data of the Crab Nebula are marked as circle to distinguish them from the others, (bottom) The ratios between the data 
to the best-fit models are shown. 
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The spectrum outside i?3% is fitted with a phenomeno- 
logical model of diskbb + powerlaw + gaussian, keep- 
ing the center energy of the gaussian fixed at 6.4 IceV. The 
resultant parameters are Tin — 0.22 keV, normalization of 
diskbb of 4.8 xlO^, F = 2.0, normalization of powerlaw 
of 3.6 photons keV~^ cm~^ second"^, and an equivalent 
width of 264 eV. The spectral ratio of the whole region 
(blue) to the model shows typical pileup features similar 
to the ratio in figure 2; i.e., the moderate hardening, the 
apparent decrease in flux, and the spectral steeping at the 
highest end around ^7 keV. On the other hand, the spec- 
tral ratio taken outside Ri% (red) agrees with that outside 
i?30/„ (black) within ^ 5% above ^ 2 keV. Thus, the spec- 
trum outside i?3% is almost free from piled-up events. 
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■*-*■***-■"* *** **>V.t . -lit it i 



12 5 10 

Energy (keV) 



Fig. 7. (top) The detector-response-convolved spectra of the 
XISO of Cyg X-1 acquired on June 6 in 2009, accumulated 
from the region outside R3% (46 pixel) and Ri% (88 pixel), 
presented in black and red, respectively. As a references, the 
spectra taken from the whole region is also shown in blue, 
(middle) The same as top ones except for the detector-re- 
sponse-deconvolved spectra shown in a uFu form, (bottom) 
The ratio to the model of diskbb + powerlaw + gaussian, 
determined with the black spectrum. 



5.3. Other systematic uncertainties 

For mitigating pileup effect, the larger mask radius is 
the better. However, at a large radius, we encounter dif- 
ferent systematic uncertainties. In order to investigate 
these, we used the dimmest (OBSID=402072020) among 
the 25 Cyg X-1 observations (see Yamada 2011 for de- 
tails). For the source, and Ri% of XISO and XISl 
are 8 pixel and 35 pixel, and 10 pixel and 38 pixel, re- 
spectively. To illustrate systematic errors outside of 
we added i?o.i%, which corresponds to ^ 100 pixel, to as 
the mask radius and showed in figure 8. The spectra out- 
side 8 pixel (i?3%) are fitted with a model of diskbb + 
powerlaw + gaussian, which is compared to the other 
spectra. 

Although all the spectra above ~ 2 keV almost agree 



with each other within ^ 5-10 %, while the spectra taken 
outside of i?o.i% (green) deviate by ~ 20% around 0.7-2 
keV only for the XISO. Thus, it is likely that the deviation 
is caused by incomplete modeling of of the spatial distri- 
bution of the contamination accumulating on the CCD 
surface, because of thicker contamination in XISO than 
XIS3''. This illustrates that a larger mask radius for a 
smaller pileup fraction is not always better in spectral 
analysis. We thus recommend to use either i?3% or Ri% 
as a practical solution. 

6. Summary 

We have comprehensively studied pileup phenomena of 
the XIS; by analyzing the actually observed data for sev- 
eral bright point sources. Using the surface brightness, the 
pileup fraction as the primary measure, we characterized 
grade-branching ratios of split events and detached events. 
We examined the effects on spectral analysis when exclud- 
ing the central regions within i?3% and Our results 
are summarized as follows: 

• Pileup causes spectral hardening in FI and BI CCDs, 
as well as changes in the sharp spectral features 
such as atomic edges of instrumental response. The 
changes in hardness is not dependent on intrinsic 
spectral hardness. 

• Pileup increases the ratio of split events by ^ 20% 
from a pileup fraction of ~ 1% to ~ 10% in both 
FI and BI CCDs. On the other hand, the ratio 
of the detached events increases proportionally with 
the pileup fraction by several orders of magnitude 
from ^0.1 and ~ 10 % of pileup fraction. 

• When surface brightness, spectral hardness, and 
split event ratio begin to increase, the pileup frac- 
tion starts to exceed ~ 1%, which illustrate that the 
pileup fraction can be a good measure for estimating 
pileup. 

• The relation between the measured surface bright- 
ness and the radius corresponding to pileup fraction 
of 3% and 1%, is useful to estimate the extent of 
pileup in a practical manner. 

• Based on the observation of Cygnus X-1, the differ- 
ences between the spectra accumulated outside i?3% 
and Ri% are less than 5% in the 2-10 keV and ~ 
20 % in 0.5-2 keV. The latter might be caused by 
uncertainties of contamination modeling in the XIS. 

Thus, we have phcnomenologically compiled the pileup 
effects specific to the XIS + XRT systems and found a 
practical way of estimating pileup effects, which would 
be useful for those who analyze the bright sources taken 
with Suzaku, and even with similar CCD sensors to the 
XIS-f-XRT system, such as Soft X-ray Imager on board 
ASTRO-H (Tsunemi et al. 2010). 

The authors would like to express their thanks to 



^ See Figure 7. 17 of Suzaku technical description for A07. 
http: / /www. astro. isas.ac.jp/suzaku/proposal/ao7/ 
suzaku_td_ao7.pdf 
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Fig. 8. The same as figure 7, but for the observation of Cyg X-1 taken in May 17th on 2007, and 
for (a) XISO and (b) XIS3 data. The spectra taken outside ^?o.i% (~100 pixel) are added in green. 



